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Abstract We study the localized direct CP violation 
in the hadronic decays —>■ p°(a;)7r^ —>■ 

including the effect caused by an interesting mecha¬ 
nism involving the charge symmetry violating mixing 
between and w. We calculate the localized integrated 
direct CP violation when the low invariant mass of 
TT+TT” [m{'K^TT~)iow\ is near p°(770). For five models 
of form factors investigated, we find that the localized 
integrated direct CP violation varies from -0.0752 to 
-0.0290 in the ranges of parameters in our model 
when 0.750 < m{‘n^TT~)iow < 0.800GeV. This result, 
especially the sign, agrees with the experimental data 
and is independent of form factor models. The new 
experimental data shows that the signs of the localized 
integrated CP asymmetries in the regions 0.470 < 
m{'K^-K~)iow < 0.770GeV and 0.770 < m['K^-K~)iow < 
0.920 GeV are positive and negative, respectively. We 
hnd that p-w mixing makes the localized integrated 
CP asymmetry move toward the negative direction, 
and therefore contributes to the sign change in those 
two regions. This behavior is also independent of 
form factor models. We also calculate the localized 
integrated direct CP violating asymmetries in the 
regions 0.470 < m{TT~^Tr~)iow < 0.770GeV, 0.770 < 
m(7r+7r“);ou) < 0.920 GeV and the whole region 0.470 < 
m{'K^-K~)iow < 0.920 GeV and find that they agree with 
the experimental data in some models of form factors. 
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1 Introduction 

GP violation is one of the most fundamental and 
important properties of weak interactions. Even though 
it has been known since 1964 [I], we still do not know 
the source of GP violation completely. In the Standard 
Model, CP violation originates from the weak phase 
in the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
m- Besides the weak phase, a large strong phase 
is also needed for direct CP violation to occur in 
decay processes. Usually, this large phase is provided 
by QCD loop corrections and some phenomenological 
mechanisms. In the past few years, numerous theo¬ 
retical studies have been conducted on CP violation. 
However, we need a lot of data to test these approaches 
because there are many theoretical uncertainties such 
as CKM matrix elements, hadronic matrix elements, 
and nonfactorizable effects. These uncertainties would 
be reduced by the increase of experimental data in the 
future and the improvement of theoretical methods. 

Recently, the LHCb Collaboration focused on three- 
body final states in the decays of B and D mesons 
and a novel strategy to probe CP asymmetry in their 
Dalitz plots [115116] . The local asymmetries in specific 
regions of the phase space of charmless three-body 
decays of bottom mesons, such as B^ 7r=*=7r+7r“ and 
B^ —?► A^7r+7r“, were measured. It was shown that 
the local asymmetry distributions in the Dalitz plots 
reveal rich structures and are not uniform |1[51I5] . These 
intriguing discoveries offer opportunities to search for 
different sources of CP violation, through the study of 
the signatures of these sources in certain phase spaces 
of the Dalitz plots. In fact, several theoretical studies 
have been made to explain these distributions, such as 
the interference between intermediate states [7115115] and 
final-state rescatterings [iiiniiiiiiii]. One can confirm 







2 


that these complex structures originate from more than 
one source [B]. 

Charge symmetry is broken at the most fundamen¬ 
tal level in strong interaction physics through the small 
mass difference between up and down quarks in the 
QCD Lagrangian. As a consequence the physical and 
w mesons are not eigenstates of isospin but, for example, 
the physical contains a small admixture of an / = 0 
state m- In previous works, this phenomenon, known 
as p-Lu mixing, was considered to obtain a large strong 
phase in B decays [HnsiiiB] and it was found that 
such a mixing can lead to a peak of CP violation when 
the invariant mass of 7 r"'' 7 r“ is near to. The differential 
CP asymmetry was studied in the decays ^ 

p°(aj) 7 r=*= — >• 7 r+ 7 r“ 7 r=*= before HZllllllS]. p-oj mixing 
provides additional complex terms to the amplitudes 
PP] and the strong phase passes through 90 ° at 
the w resonance [niiainiiiiiiiiii]. In our previous 
work, we found that it is more useful to investigate 
the localized integrated CP asymmetry and studied 
the localized integrated CP asymmetry in hadronic r 
decays with this mechanism [20]. The newest LHCb 
experiments showed that the resonances are /o°(770) 
when TO^(7r+7r“)iou) < IGeV^/c^ for B^ —>■ 
decays [m( 7 r“'' 7 r“);ou, is the low invariance mass of 
7 r+ 7 r“] [Bj. It was claimed that the CP asymmetry in 
these decays changes sign around the p°{770) peak 
of |B], which contains the significant 

region of p-uj mixing. In the present paper, we aim 
at studying the localized integrated CP asymmetry 
in B^ — >■ TT+TT^TT^ decays involving p-ui mixing and 
comparing it with the results obtained from the LHCb 
Collaboration. 

In this paper, we will investigate the localized inte¬ 
grated CP violation with five phenomenological models 
of weak form factors for B^ —>■ p^{uj)n^ —>■ tt^tt+tt^ 
with and without p-uj mixing. We will determine the 
allowed range of W which is the effective color number 
describing nonfactorizable contributions. The model 
dependence of our results will be discussed in detail. 
We will see that for five models of form factors, our 
result for the localized integrated CP asymmetry varies 
from -0.0752 (-0.0626) to -0.0403 (-0.0290) in the ranges 
of Nc, 2.07(2.09) < Nc < 4.54(4.65) corresponding to 
q^lml = 0.3(0.5) {q is the typical momentum transfer 
of the gluon or photon in the penguin diagrams) 
when 0.750 < m(7r+7r“);ou) < 0.800 GeV. From Fig. 4 
of Ref. |B], we can see the localized integrated CP 
asymmetries have different signs in cos 0 > 0 and 
COS0 < 0 regions when 0.750 < m(7r+7r“);ou) < 
0.800 GeV {9 is the angle between the momenta of 
the unpaired tt and the p^ decay product with 
the same-sign charge). If one adds the events in these 


two experimental regions together, the total localized 
integrated CP asymmetry will be —0.0294 ± 0.0285. 
Our results agree with this experimental data. The 
experimental values of the localized integrated CP 
asymmetries in the regions 0.470 < m{'!T'^Tr~)ioTju < 
0.770GeV and 0.770 < m( 7 r+ 7 r-);o^ < 0.920GeV are 
0.0508 ± 0.0171 and —0.0256 ± 0.0202, respectively, 
with opposite signs |B]. We will find that p-ui mixing 
can make the localized integrated CP asymmetry move 
toward the negative direction, and therefore will con¬ 
tribute to the sign change from the region 0.470 < 
m(7r+7r“)iou) < 0.770GeV to 0.770 < m{'K^'n~)iow < 
0.920 GeV. This behavior is independent of form factor 
models. Furthermore, we will find that our results 
in these two regions and the whole region of them 
are consistent with the experimental data for several 
models of form factors. 

The remainder of this paper is organized as follows. 
In Sect. 2, we present the form of the effective Hamilto¬ 
nian and the values of Wilson coefficients and give the 
formalism for the CP violating asymmetry in B^ —>■ 
P^(uj)tt^ —>■ TT+TT^TT^. Then we show numerical results 
with several models of form factors in this section. In 
Sect. 3, we calculate branching ratios for 5+ —>■ p° 7 r+ 
and B^ —>■ p+ 7 r“ and present numerical results for the 
range of Nc allowed by the experimental data. In the 
last section, we give some discussions and summarize 
our results. 


2 CP violation in —>■ p°(ci;)7r=*= —>■ tt+tt 

The amplitude of a decay process described by some 
amplitudes may have CP-even and -odd relative phases. 
Within the Standard Model, the CP-odd relative phase 
is always a weak phase difference which is directly 
determined by the CKM matrix. On the other hand, 
CP-even phases, which are called strong phases, usu¬ 
ally originates from nonperturbative effects of strong 
interactions and are hard to handle. We consider a 
B meson weak decay process, B —>■ Mi M 2 M 3 , where 
Mi(f = 1, 2,3) is a light pseudoscalar meson. For a weak 
decay process of a heavy meson, a typical form of the 
decay amplitude A and its CP conjugate one A are 

A = girie“^i -f (1) 

A = -b 52^2e“^^ (2) 

where gi and 52 represent CP-odd complex terms which 
involve CKM matrix elements, and r 2 e'^'^ terms 

are even under the CP transformation. Then one has 

|Ap - |Ap = 4rir2 '3m{glg2) sin((()i - (f>2) 

= 4 rir 2 |gi||g 2 | sin[Arg( 52 /ffi)] sm{(j3i - (j) 2 ), (3) 
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from which, we can see explicitly that both the CP- 
odd phase difference Arg(( 72 /ffi) and the CP-even phase 
difference (pi — (j )2 are needed to produce CP violation. 


2.1 The effective Hamiltonian 

In order to calculate the direct CP violating asymmetry 
in hadronic decays, one can use the following effective 
weak Hamiltonian, based on the Operator Product 
Expansion 


„ 10 

n = ^ [Vui,v:,{ci0^+c20^)-vM E 

(4) 


where Vub, Vud, Vtb and Vtd are CKM matrix elements, 
and Ci{i = 1, 2,.., 10) are the Wilson coefficients, which 
are calculable in the renormalization group improved 
perturbation theory and are scale dependent. In the 
present case, we work with the renormalization scheme 
independent Wilson coefficients and use the values of 
the Wilson coefficients at the renormalization scale n ~ 
mb- The operators Oi have the following forms: 

O" = - 75)U/3U/37'^(1 - 75)&a, 

O 2 = ^7^(1 - 75)'UU7'"(1 - 75 )^ 1 , 

O 3 = ^7^(1 - - 75)<?', 

q' 

O 4 =rfa7M(l-75)^/3E^/3T''"(l“T'5)'Ia> 
q' 

O 5 = ^7^(1 - 75)&E9'T'''(4 + T'5)g', 
q' 

Oe =da7M(l-75)fe/3E^/3"l''(4+T'5)9a> 

q' 

O 7 = + J 5 )q', 

q' 

Os = - "ib)bl3 E + l5)q'a, 

q' 

O9 = ^d 7 A«(l- 75 )&E® 9 '? 7 ^(l- 75 )g', 
q' 

OlO = - lb)bl3 E - lb)q'a., (5) 

9' 

where a and /3 are color indices, and q' = it, c? or 
s quarks. In Eq. (O, 0“ and O 2 are the tree level 
operators, O 3 — Oq are QCD penguin operators, and 
O 7 — Oio arise from electroweak penguin diagrams. 

The Wilson coefficients, c,, which are known to the 
next-to-leading logarithmic order, take the following 
values HZUIH]: 


Cl = -0.3125, C 2 = 1.1502, 

C 3 = 0.0174, C 4 = -0.0373, 

C 5 = 0.0104, C 6 = -0.0459, 

C 7 = -1.050 X 10-^ C8 = 3.839 x 10-^ 
eg = -0.0101, CIO = 1-959 X 10"^ 


( 6 ) 


at the scale /r = mf, = 5 GeV. 

To be consistent, the matrix elements of the oper¬ 
ators Oi should also be renormalized to the one-loop 
order. This results in the effective Wilson coefficients, 
c', which satisfy the constraint 


Ci{mb){0^{mb)) = c-(Oi)*’'®®, 


(7) 


where (Oi)*’'®® is the matrix element at the tree level, 
which will be evaluated in the factorization approach. 
Erom Eq. 0 , the relations between c' and Ci are [27l 

[2H1 


Cl — Cl, C2 

C 3 = C 3 - -Ps/3, 
4 = C 5 - -Ps/3, 

C7 = C7 -|- He, 

Cg = Cg -\- Pei 

where 


C2, 

C4 = C4 -|- Hjj, 
4 = Co -f Hs, 
Cg = Cg, 

Cio — ClO, 


( 8 ) 


Ps = (as/87r)c2[10/9-bG(mc,M, g^)], 

Pe — (o^em /97r)(3ci -I- C2)[10/9 -I- G{mc, fi, g^)], 


with 

G(mc,fJ,,q ) = 4 / dxx(x — l)ln---. 

Jo M 

G{mc, fJ.,q^) has the following explicit expression 



Based on simple arguments at the quark level, the value 
of g^ is chosen to be in the range 0.3 < g^/m^ < 0.5 
Erom Eqs. ( 8 ) and (9) we can obtain numerical 
values for c'. 
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When Iml = 0.3, 

c'l = -0.3125, 4 = 1.1502, 

4 = 2.433 X 10"^ + 1.543 x 10"^i, 

4 = -5.808 X 10-2 - 4.628 x W^i, 

4 = 1.733 X 10-2 + 1.543 X-3 i, 

4 = -6.668 X 10-2 - 4.628 x-^ i, 

4 = -1.435 X 10-'‘ - 2.963 x-^ i, 

4 = 3.839 X 10-4 
4 = -1.023 X 10-2 - 2.963 x-® i, 

4o = 1.959 X 10-4 (10) 

and when /ml = 0.5, one has 

4 = -0.3125, 4 = 1.1502, 

4 = 2.120 X 10-2 + 2.174 X lO-^i, 

4 = -4.869 X 10-2 - 1.552 x 10-^i, 

4 = 1.420 X 10-2 + 5.174 X-3 i, 

4 = -5.729 X 10-2 - 1.552 x -2 i, 

4 = -8.340 X 10-^ - 9.938 x-® i, 

4 = 3.839 X 10-4 
4 = -1.017 X 10-2 - 9.938 X-® i, 

4o = 1.959 X 10-4 (11) 

where we have taken as{mz) = 0.112, aem{mb) = 
1/132.2, mb = 5GeV, and = 1.35 GeV. 


2.2 Formalism 


For the B —>■ p^n process, the amplitude can be 
written as = apB ■ e*(A), where e is the 

polarization vector of and A is its polarization, ps 
is the momenta of B~ meson, and a is independent 
of A. The amplitude for —>• tt+tt- is = 

gpe{X){pi — P 2 ), where pi and p 2 are the momenta of 
7r+ and tt~ produced by respectively, and Pp is the 
effective coupling constant for —>■ tt+tt-. Then, for 
the sequential decay B~ —>• p^Tr~ —^ tt+tt-tt-, the 
amplitude is [511^] : 


A = 


_ p 


ap'g- 


9p^ p 


9p{pi-V2Y 


9pp 


{pi +P2)p[Pl +P2)i^ 


[Pl - P 2 T 


_ 9p ^B+^p°Tr+ / /N 

“ 7-1--1^ - S J 


Pb ■ e 


= {E-s')-M, 


( 12 ) 


where a/s' is the high invariance mass of the 7r"''7r 
pair, r = \{aYax + s'min) with s'^ax and s4„ being 
the maximum and minimum values of s' for a fixed 


s, respectively, and 4® is the low invariant mass of 
the Tr+Tf- pair [TO(7r''"7r-);ou,] and sy is from the 
inverse propagator of the vector meson V, sy = 
s — my + irnyBy. The factor E — s' is equal to 
—2 |p2||p 3| COS0 in Ref. [30] , which accounts for angular 
momentum conservation for the spin-1 resonance. M. 
will be calculated in the following. According to the 
effective Hamiltonian, A can be divided into two parts 

m- 

A = 'k~\H'^\B~) + {'k''^tt~'k~\H^\B~)^ (13) 

where and are the Hamiltonians for the tree 
and penguin operators, respectively. 

In order to obtain a large signal for direct CP 
violation, we need to appeal to some phenomenological 
mechanisms, p-uj mixing has the dual advantages that 
the strong phase difference is large (passes through 90 ° 
at the w resonance) and well known [TSlll6j . With this 
mechanism, to the first order in isospin violation, the 
amplitude for B~ —>■ p^{u!)Tr~ —>■ tt+tt-tt- takes the 
following form at a value of ^/s close to the uj resonance 
mass 

{TT+n-Tr-\H^\B-) = (A - s') 

xf-^flp^t^ + ^tp) , (14) 

\SpSt^ Sp J 

(7r+7r-7r-|i7^|H-) = (A - s') 

X (-^^flpujPuj + —Pp) , (15) 

\SpSui Sp J 

where ty {V = or w) is the tree amplitude 
and py is the penguin amplitude for producing an 
intermediate vector meson V, Hp^ is the effective p- 
bj mixing amplitude. From Eqs. (14) and (15), we note 
that p-Lo mixing provides an additional complex term 
for the tree and penguin amplitudes (the first term 
in each equation), respectively. These complex terms 
will enlarge the CP-even phase and lead to a peak 
of CP asymmetry as mentioned before. We will show 
the difference between the CP asymmetries with and 
without p-uj mixing later. Here, we assume that the 
—> 7r''’7r TT^ process is dominated by the resonance 
in certain region of its Dalitz plot. 

We stress that the direct coupling w —)► tt+tt- 
is effectively absorbed into TIpui EH, leading to the 
explicit s dependence of iTpoj. Making the expansion 
npu{s) = npuj{ml,) + {s-ml/)n'pYBni), the p-w mixing 
parameters were determined in the fit of Gardner and 
O’Connell [5^ : 

meUp^irnl) = - 3500 ± 300 MeV2, 

3vn.npi^imY) = — 300 ± 300 MeV2, 
n'p^{ml) = Qm± 0.04. 


( 16 ) 











5 


In practice, the effect of the derivative term is negligi¬ 
ble. 


In this work, we only consider and w resonances. 
Then, for a fixed s, the differential CP asymmetry 
parameter can be defined as 
|7W|2-|7W|2 


Acp = 


\MV + |7Wp 


(17) 


By integrating the denominator and numerator of Aqp, 
respectively, in the region 17 (si < s < S 2 , si < s' < 
S 2 ), we obtain the localized integrated CP asymmetry, 
which can be measured by experiments and takes the 
following form: 


ID f _ O' 

/;;ds/;;=ds'(i:-s02(iMp4-|M|^) 


(18) 


According to kinematics of the three body decay, IJ[= 
^Wmax + Smm)] related to s. In our calculations, 
s varies in a small region, and therefore E can be 
treated as a constant approximately [B]. Then, the 
terms ds'(i7 —s')^ are cancelled, and Afip becomes 
independent of the high invariance mass of tt+tt”. In 
practice, to be more precise, we take into account the 
s-dependence of and s'^^^ in our calculations. 

We choose s'^^^ < s' < as the integration 

interval of the high invariance mass of tt+tt” and regard 
jj>max ^ factor which is dependent on s. 


2.3 Calculational details 

With the Hamiltonian given in Eq. (4), we are ready 
to evaluate the matrix elements for B~ —>■ p^{uj)n~. 
According to the theory of QCD factorization, the 
naive factorization approach has been shown to be 
the leading order result in the framework of QCD 
factorization when the radiative QCD corrections of 
order 0{as{mb)) and the 0{l/mb) corrections in the 
heavy quark effective theory are neglected miMiiisi 
I26| . Since the b quark is very heavy and B meson decays 
are very energetic, so the quark-antiquark pair in a 
meson in the final state moves very fast away from the 
weak interaction point. The hadronization of the quark- 
antiquark pair occurs far away from the remaining 
quarks. Then the meson can be factorized out and the 
interaction between the quark pair in the meson and 
the remaining quark is tiny PIIM] . The deviation of 
the value of from the color number, 3, measures the 
nonfactorizable effects in the naive factorization scheme 
[min]. In the factorization approximation, either p^{uj) 
or 7 r“ is generated by one current which has appropriate 
quantum numbers in the Hamiltonian. For this decay 
process, two kinds of matrix element products are 


involved after factorization: (p°(a;)|J''| 0 )( 7 r“| 
and (7r“|J^|0)(p°(a;)|We will calculate them 
in some phenomenological quark models. 

The matrix elements for H P and B ^ V (where 
P and V denote pseudoscalar and vector mesons, 
respectively) can be decomposed as [35] 


{P\J^\B) = (^PB +pp 


uid — m‘ 




/ a 


fc2 


-k^F^ik^) 


{V\J^B) = 


--- {k'^) 

rriB -I- my 

-bi|e* (ms + my)Af ^(/c^) 
e* • k 


rriB -b my 
e*-k 


A* 

BV li„2 


k^ 


2 my • kpA^'^ [k 


+i^^2my • kpAQ^ik"^), 


(19) 


{pb + Pv)tJ.Ai^ {k'^) 


( 20 ) 


where Jp is the weak current [J^ = — 75)6 or 

— k = Pb —pp{v)^ a-nd is the polarization 
vector of V. The form factors included in our calcula¬ 
tions satisfy Pf^(O) = Po^^(O), Af^(O) = A^'^(O), 
and A^^(k"^) = [(ms -b my)/2my]Af ^(fc^) — [(ms — 
my)/2my]Af'^(fc^). We define the notation X for 
matrix elements. For example, ^ Ms defined as 
( 7 r“|J 7 ''(l — 75 )it1 0 )(p°|? 27 ^ (1 — 75)61 P“). These matrix 
elements can be parameterized as the products of decay 
constants and form factors. Therefore, the factorized 
terms have the following expressions [36] : 


X^^FV) ^ - 75)g3|0)(P|gi7Ml - l5)b\B) 

= ‘2fvmvF^^{my){e* ■ pb), ( 21 ) 

xiBV,p) ^ (p|^ 27M1 - 75)g3|0)(F|gi7Ml - l5)b\B) 

= 2/pmyA^^(mp)(e*-ps), (22) 


where fy and fp are the decay constants of vector and 
pseudoscalar mesons, respectively. Using the decompo¬ 
sition in Eqs. (19)-(22), one has 


tp — ^ubVy^ii 

M I 


e* - PB 

j^{B-ui,n-) 

e-*-PB 




M) 


0,2' 


02' 


- PB 

e* - PB 


(23) 

(24) 


where all the Oi are built up from the effective Wilson 
coefficients c'’s, and take the form Oi = c' -bc'_|_i/A^c for 
odd i and Oi = c'-bcM^/A^c for even i. It is noted that in 
the factorization approach Nc includes nonfactorizable 
contributions effectively, and the value of Nc should be 
determined by experiments since we cannot evaluate 
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nonfactorizable contributions. In the same way, we 
obtain the penguin operator contributions: 


Pp = -VM 

+ 04 + Oio — 2(05 + Og) 


3 3 1 

— 04 + + 2®9 + 2 *^^° 


X{B TT ,p°) 

e* -PB 


{md + TO„)(mb + rriu)- 


X{B p”,7r ) 

e* • PB 


(25) 


Pu = -VtbVt 


td ' 


203 + 04 + 205 + 2(^7 + O9 ~ Olo) 

X{B-Tr-,u:) 


+ 


e* -pB 

O 4 + Oio — 2 (o 6 + Og) 
X(B-u),Tr-) I 


{md + mu){mb + m„)- 


(26) 


We adopt the same decay constants and form fac¬ 
tors for the matrix elements producing and w 
mesons. Then we have '> = ^ = 

2f^mpAo{ml){e* -ps) and = 

\^fpmpFi{ml){e* ■ pb), where 
(p°(a;)| J'^IO) = l/V2fpmpe*P and (7r“| J^|0) = if^Pp- 


2.4 Numerical results 

In our numerical calculations we have several param¬ 
eters: q^, Nc, and the CKM matrix elements in the 
Wolfenstein parametrization. As mentioned in Sect. 2.1, 
the value of is conventionally chosen to be in the 
range 0.3 < q^/ml < 0.5. The CKM matrix, which 
should be determined from the experimental data, 
has the following form in terms of the Wolfenstein 
parameters. A, A, p, p m- 

/ I — A AA^(p—i77)\ 

-A 1 - iA2 AA^ , (27) 

\ AA^(1 — p — ijy) —AA^ 1 / 

where O(A^) corrections are neglected. Since A and A 
are well determined and the uncertainties due to the 
CKM matrix elements are mostly from p and 77 , we 
take the central values of A(= 0.225) and A(= 0.814) 
in the following. The ranges for p and p are 

p = 0.117± 0.021, 77 = 0.353 ±0.013, (28) 


with 



The form factors Fi{m?p) and Ao(777,^) depend on the 
inner structure of the hadrons. Under the nearest 
pole dominance assumption, we take the following 
dependence of the form factors: for Model 1(2) [3^IM1 
[^I4ni^l4^ : 

Fiik^) = Mk^) = (29) 


where hi = 0.25(0.292), Hao = 0.30(0.366), 

777-1 = 5.32 GeV, and toaq = 5.27GeV; for Model 3(4) 


Fi{k^) 


(1 


hi 

/ 


Ao(fc2) 


It-Ao 




(1-^) 

\ mi._ / 


(30) 


where hi = 0.25(0.292), hAg = 0.30(0.366), 7774 = 
5.32 GeV, and mAo = 5.27 GeV; for Model 5 [iHlHH] : 


Fi{e) 

Ao(fc2) 


mf \mf ) 

hAo 




(31) 


where hi = 0.261, /iaq = 0.302, xi = 2.03, yi = 1.29, 
xq = —1.49, 7/0 = 6.61, 7771 = 5.32 GeV, and mAp = 
5.27 GeV. The decay constants used in our calculations 
are fp = 216 MeV and U = 132 MeV @3]. 

In the numerical calculations, it is found that 
for a fixed there is a maximum value for the 
differential CP violating parameter, when the low 
invariant mass of the tt+tt” pair is in the vicinity of 
the uj resonance, 0.780 — 0.785 GeV. Five models with 
different form factors were investigated to study the 
model dependence of Acp in Ref. m- To be more 
specific, in Figs, la, b we display the results for the form 
factors in Model 1. These results show explicitly the 
dependence of the GP violating asymmetry on 9 ^/ 777 ^, 
CKM matrix elements and the effective parameter Nc- 
The dependence on Nc comes form the fact that W is 
related to the hadronization effects, and consequently, 
we cannot exactly determine W in our calculations. 
Therefore, we treat Nc as a free effective parameter and 
take it in the range 2.07(2.09) < Nc < 4.54(4.65) when 
q^/m/ = 0.3(0.5) for reasons which will be explained 
later (Sect. 3). 

Then we calculate the localized integrated CP asym¬ 
metries. According to Eq. (18), we integrate Acp over 
the low invariant mass of 7r"''7r“ (i/s) and obtain the 
localized integrated asymmetries A^p. Considering the 
significant region of p-uj mixing, we choose the integra¬ 
tion interval of ^/s to be from 0.750 to 0.800 GeV. In 
order to compare with the newest result of the LHCb 
experiments, we also calculate A^p when ^/s is in the 
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Table 1 The localized integrated asymmetries A^p when = 0.3(0.5). The first and second lines of each model 

corresponds to A^p with and without p-uj mixing in the region 0.750 < y/s < 0.800 GeV, respectively. The third and forth lines 
of each model correspond to the low-mass region (0.470 < y/s < 0.770 GeV) and the high-mass region (0.770 < y/s < 0.920 GeV) 
near the resonance mass, respectively. The fifth and sixth lines of each model corresponds to A^p with and without p-uj mixing 
in the region 0.470 < y/s < 0.920 GeV, respectively. The second and third columns correspond to = 2.07(2.09). The fourth 
and fifth columns correspond to Nc = 4.54(4.65). The second and fourth columns correspond to lower limiting values of the 
GKM matrix elements. The third and fifth columns correspond to upper limiting values of the GKM matrix elements. 


N, 

2.07(2.09) 

4.54(4.65) 

P,V 

min 

max 

min 

max 

Model 1 





0.750 < < 0.800 GeV 

0.470 < Vs < 0.770 GeV 
0.770 < Vs < 0.920 GeV 
0.470 < Vs < 0.920 GeV 

Model 2 

-0.0647(-0.0520) 

0.0054(0.0204) 

-0.0024(0.0061) 

-0.0352(-0.0128) 

-0.0167^0.0021) 

0.0054(0.0204) 

-0.0736(-0.0591) 

0.0062(0.0232) 

-0.0027(0.0070) 

-0.0400(-0.0145) 

-0.0189^0.0023) 

0.0062(0.0232) 

-0.0455(-0.0300) 

0.0079(0.0281) 

0.0017(0.0278) 

-0.0239(0.0014) 

-0.0093(0.0096) 

0.0079(0.0281) 

-0.0517(-0.0341) 

0.0090(0.0319) 

0.0019(0.0290) 

-0.0272(0.0016) 

-0.0106(0.0109) 

0.0090(0.0319) 

0.750 < Vs < 0.800 GeV 

0.470 < Vs < 0.770 GeV 
0.770 < Vs < 0.920 GeV 
0.470 < ^ < 0.920 GeV 

Model 3 

-0.0661(-0.0516) 

0.0062(0.0234) 

-0.0019(0.0086) 

-0.0359(-0.0112) 

-0.0166(0.00004) 

0.0062(0.0234) 

-0.0752(-0.0587) 

0.0071(0.0267) 

-0.0021(0.0098) 

-0.0408(-0.0127) 

-0.0189(0.00004) 

0.0071(0.0267) 

-0.0469(-0.0290) 

0.0092(0.0327) 

0.0026(0.0306) 

-0.0247(0.0038) 

-0.0091(0.0129) 

0.0092(0.0327) 

-0.0533(-0.0330) 

0.0105(0.0371) 

0.0030(0.0353) 

-0.0280(0.0043) 

-0.0103(0.0146) 

0.0105(0.0371) 

0.750 < Vs < 0.800 GeV 

0.470 < Vs < 0.770 GeV 
0.770 < Vs < 0.920 GeV 
0.470 < Vs < 0.920 GeV 

Model 4 

-0.0647(-0.0520) 

0.0054(0.0204) 

-0.0024(0.0061) 

-0.0350(-0.0125) 

-0.0167^0.0021) 

0.0054(0.0204) 

-0.0737(-0.0592) 

0.0062(0.0232) 

-0.0028(0.0069) 

-0.0398(-0.0142) 

-0.0190^0.0023) 

0.0062(0.0232) 

-0.0455(-0.0300) 

0.0079(0.0281) 

0.0017(0.0278) 

-0.0238(0.0016) 

-0.0094(0.0097) 

0.0079(0.0281) 

-0.0515(-0.0341) 

0.0090(0.0320) 

0.0019(0.0292) 

-0.0270(0.0016) 

-0.0106(0.0110) 

0.0090(0.0320) 

0.750 < Vs < 0.800 GeV 

0.470 < Vs < 0.770 GeV 
0.770 < Vs < 0.920 GeV 
0.470 < Vs < 0.920 GeV 

Model 5 

-0.0661(-0.0517) 

0.0062(0.0235) 

-0.0019(0.0086) 

-0.0357(-0.0110) 

-0.0167(0.00005) 

0.0062(0.0235) 

-0.0752(-0.0587) 

0.0071(0.0267) 

-0.0022(0.0098) 

-0.0405(-0.0125) 

-0.0190(0.00005) 

0.0071(0.0267) 

-0.0469(-0.0290) 

0.0092(0.0327) 

0.0026(0.0306) 

-0.0245(0.0040) 

-0.0091(0.0129) 

0.0092(0.0327) 

-0.0534(-0.0330) 

0.0105(0.0372) 

0.0030(0.0353) 

-0.0278(0.0046) 

-0.0104(0.0146) 

0.0105(0.0372) 

0.750 < Vs < 0.800 GeV 

0.470 < Vs < 0.770 GeV 
0.770 < ^ < 0.920 GeV 
0.470 < Vs < 0.920 GeV 

-0.0577(-0.0550) 

0.0011(0.0041) 

-0.0055(-0.0083) 

-0.0322(-0.0216) 

-0.0170^0.0140) 

0.0010(0.0039) 

-0.0658(-0.0626) 

0.0012(0.0047) 

-0.0063(-0.0094) 

-0.0366(-0.0245) 

-0.0194(-0.0160) 

0.0012(0.0044) 

-0.0403(-0.0375) 

0.0015(0.0055) 

-0.0034(-0.0046) 

-0.0222(-0.122) 

-0.0114(-0.0079) 

0.0015(0.0052) 

-0.0459(-0.0427) 

0.0017(0.0062) 

-0.0039b0.0053) 

-0.0252(-0.0139) 

-0.0130b0.0090) 

0.0017(0.0059) 


low-mass region (0.470 < y/s < 0.770GeV), the high- 
mass region (0.770 < y/s < 0.920 GeV) and the whole 
region (0.470 < y/s < 0.920 GeV) near the resonance 
[5]. The numerical results are displayed in Table 1. We 
also display Aqp with and without p-ui mixing when 
0.750 < yi < 0.800 GeV and 0.470 < < 0.920 GeV 

in Table 1. 

Table 1 shows that the values of A^p in our 
calculations vary from -0.0752 (-0.0626) to -0.0403 (- 
0.0290) corresponding to q^/ml = 0.3(0.5), in the 
regions of the GKM matrix elements, and the 
form factors in five models when 0.750 < y/s < 
0.800GeV. From Fig. 4 of Ref. [5], we can see the 
localized integrated CP asymmetries have different 


signs in cos 0 > 0 and cos 0 < 0 regions when 
0.750 < m(7r+7r“)iou, < 0.800GeV. If one adds the 
events in these two experimental regions together, the 
total localized integrated CP asymmetry obtained from 
experiment becomes —0.0294 ± 0.0285 when 0.750 < 
m{TT^Tr~)iow < 0.800 GeV. The values in our calcu¬ 
lations agree with this experimental data. We stress 
that Aqp in our calculations is always negative in 
this integration region and its sign is independent of 
form factor models. We note that the signs of A^p 
are positive when p-uj mixing is not considered in this 
region. This indicates that p-ui mixing is vital for Aqp 
to be negative in this region. 
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Fig. 1 The differential asymmetry, Acp, in Model 1. (a) For q^fm^ = 0.3: Solid and dot lines correspond Nc = 2.78 
and = 9.68 with maximum CKM matrix elements, respectively. Dot-dashed and dashed line correspond = 2.78 and 
Nc = 9.68 with minimum CKM matrix elements, respectively; (b) Same as (a) but for lm\ = 0.5 and = 2.85 and 9.02. 


From Table 1, we can see p-Lo mixing changes 
the sign of Aqp from positive to negative. Figs. 1(a) 
and 1(b) show that the peak of the differential asym¬ 
metry Acp involving p-u mixing is on the right of 
0.770 GeV. Therefore, comparing with A^p in the range 
0.470 < yi < 0.770 GeV, the localized integrated CP 
asymmetries move towards the negative direction when 
0.770 < ^/s < 0.920 GeV due to p-ui mixing. This 
behavior is also independent of form factor models. 
In fact, in our calculations the difference between the 
localized integrated CP asymmetries in the regions 
0.470 < yi < 0.770 GeV and 0.770 < < 0.920 GeV 

varies from 0.0076 to 0.0387. We add the experimental 
events for positive and negative cos 9 in the regions 
0.470 < m(7r+7r“)jouj < 0.770GeV and 0.770 < ^/s < 
0.920 GeV and obtain A^p as 0.0508 ± 0.0171 and 
—0.0256±0.0202, respectively, from the data in TABLE 
IV of Ref. [5]. We can see they have opposite signs. 
After we combine these two regions together, Aqp in 
the whole region 0.470 < m(7r+7r“);ou, < 0.920 GeV 
is 0.0173 ± 0.0130 [6]. The values of Aqp with p-uj 
mixing shown in Table 1 differ a lot between Model 1 (or 
2,3,4) and Model 5. Aqp is also sensitive to the choice 
of jm^. When 0.470 < y/s < 0.770GeV, it appears 
that Aqp varies from 0.0278 to 0.0353 in Models 1,2,3 
and 4 for q^/ml = 0.5 and Nc = 4.65. This result is 
consistent with the experimental data. When 0.770 < 
y/s < 0.920 GeV, except for q'^fml = 0.5 and W = 4.65 
in Models 1,2,3 and 4, Aqp varies from -0.0408 to - 
0.0110. This result is consistent with the experimental 


data. In the whole region of 0.470 < i/s < 0.920 GeV, 
Aqp varies from 0.0096 to 0.0146 in Models 1,2,3 and 
4 for q^lml = 0.5 and Nc = 4.65. This result is also 
consistent with the experimental data. 

From above discussions, we can see Aqp with p- 
uj mixing, especially the signs, agree with the experi¬ 
mental result when 0.750 < y/s < 0.800 GeV. We also 
find that the localized integrated CP asymmetries move 
toward the negative direction due to p-w mixing. Aqp 
in the region 0.770 < y/s < 0.920 GeV contains the 
contribution of p-uj mixing while that in the region 
0.470 < y/s < 0.770 GeV does not. Therefore, p-w 
mixing contributes to the sign change of CP asymmetry 
around the p°(770) peak of m['K~^TT~)iow Our results in 
these two regions and the whole region are consistent 
with the experimental data for several choices of q^/m? 
and models of form factors. One should take the effect of 
p-uj mixing into account in order to answer the question 
why the sign of CP asymmetry changes around the 
p°(770) peak of m(7r+7r“)jouj. 





9 


14 

13 

12 

11 

10 

o 

Ci. 

T 9 

+ 

S. 

ct1 8 
HI 

7 

6 

5 

4 


X 10" 




max. exp. value 


min. exp. value 


5 

N, 


Fig. 2 The branching ratio for B+ —^ p®7r+ in Model 2. Solid (dashed) line stands for = 0.3 and maximum (miximum) 

CKM matrix elements. Dot (dot-dashed) line stands for = 0.5 and maximum (miximum) CKM matrix elements. 


3 Extraction of TVc from data of branching 
ratios 


3.1 Formalism 


As mentioned before, includes nonfactorizable ef¬ 
fects which are difficult to deal with at present. There¬ 
fore, we treat Nc as an effective parameter to be deter¬ 
mined by the experimental data. With the factorized 
decay amplitudes, we can calculate the decay rates by 
using the following expression [36]: 


r{B VP) 


IppP A{B ^ VP) 2 
Sttuiy e* ■ Pb ’ 


(32) 


where 


^ y/[m% - {my + - {my - mp)^] 

2mB 

is the c.m. momentum of the decay particles, and 
A{B —>■ VP) is the decay amplitude. In our case, to 
be consistent, we should also take into account the p-uj 
mixing contribution when we calculate the branching 
ratio since we are working to the first order of isospin 
violation. Explicitly, for 5+ —we obtain 


BR{B+ p°TT+) 

Gl\Pp\^ 


[VudV:,AG - VuV;M/{e* -pb) 


l6nmprB+ 

+ [VudV:,Al - Vt,V,lAP]/{e*-B) 


X 


(Sp 


Ppuj 

ml) + imu^Pcj 


(34) 


where the tree and penguin amplitudes are 

3 3 1 




= ( - 04 + -07 + -09 + 

+ 04 + Oio — 2(05 + Og) 


2 


(m„ + md){mb + to„)- 


Af , — 2 o 3 + 04 + 2 o 5 + 2(^7 + 09 — Olo) 


O 4 + Olo — 2 (o 6 + Og) 

+ ■7r + 


X^B+-rr+,uj) 


{mu + md){mb + to„)- 
XX^B+U^,7r+)^ (35) 


For ^ p'^TT , we obtain 


BR{B° p+TT-) = 


g||PpP 

Wirm'^Pso 
r* A P 


x\{VubV:dAU - VtbV;dA^^)/iP-pB)\^, (36) 


where 


A^+ = (04 + Oio)x(^°^'’"^), (37) 

where ^ = 2/^TOpAo(mJ)(e* ■ pb) and 

^(B%-,P+) = ^fpmpF^{ml){e* ■ pb). 
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Table 2 Summary of the ranges of Nc which are determined from the experimental data for various models and input 
parameters [numbers outside (inside) brackets are for = 0.3(0.5)]. The notation (number, number) means the lower and 

upper limits for Nc- (-, -) means that there is no range of Nc which is consistent with the experimental data. 



B+ 

B° 

Model 1 

pmax 5 pmax 

(2.70, 4.09) [(2.76, 4.22)] 

-) [(-. -)] 

Pmin 5 Pmin 

(2.14, 3.05) [(2.18, 3.13)] 

(-. -) [(-. -)] 

Model 2 

pmax 5 pmax 

[(-.-)] 

[(-.-)] 

Pmin ! pmin 

(-r) [(-.-)] 

(-r) [(-,-)] 

Model 3 

Pmax 5 Pmax 

(2.24, 3.23) [(2.28, 3.32)] 

(0.05, 0.07) [(0.10, 0.14)] 

Pmin ! Pmin 

(2.85, 4.40) [(2.92, 4.55)] 

(0.06, 0.13) [(0.12, 0.15)] 

Model 4 

Pmax iPmax 

(-r) [(-.-)] 

(-.-) [(-,-)] 

Pmin 5 Pmin 

(-r) [(-.-)] 

(-.-) [(-.-)] 

Model 5 

pmax 5 Pmax 

(2.71,4.54) [(2.76,4.65)] 

(0.11,0.12) [(0.11,0.12)] 

Pmin 5 pmin 

(2.07,3.13) [(2.09,3.18)] 

(0.10,0.11) [(0.10,0.11)] 





3.2 Numerical results 

The latest experimental data of branching ratios of 
/9°7r+ and —>■ from Particle Data Group 

(PDG) are [37] 

BR{B+ p°7r+) = (8.3 ± 1.2) x 10"®, 

BR{B° p+TT-) = (2.3 ± 0.23) x 10"^ 

We can determine the range of W by comparing the 
theoretical values of the branching ratios with those of 
two-body decay channels. We calculate the branching 
ratios with the formula given in Eqs. (32), (34) and 
(36) in five models for the weak form factors which 
are mentioned in the previous subsection. In Fig. 2, 
we show the results for —>■ in Model 1 as 

an example. The numerical results are sensitive to 
uncertainties coming from the experimental data. In 
addition, the branching ratio also depends on the CKM 
matrix elements which are parameterized by A, A, p, 
and rj. In the allowed ranges for the parameters p and 
77 , we obtain the range of W- We summarize the allowed 
range of W in Table 2. It is found that the experimental 
data constrain the value of W into two regions. We 
note that if W approaches to zero, nonperturbative 
effects would be very large. Considering this, we drop 
the range Nc < 1. Therefore, W could be in the range 
2.07(2.09) < Nc < 4.54(4.65) for q^/ml = 0.3(0.5). 
These values have been used in Sect. 2. 


4 Conclusion and discussion 

The first aim of the present work is to study the 
localized integrated CP asymmetry for the decays 


B^ — ^ p^{uj)n^ — >■ TT+TT TT^ with the inclusion of p- 
uj mixing. The second aim is to study the sign change 
caused by p-w mixing. 

In the calculation of CP violating asymmetry pa¬ 
rameters, we need the Wilson coefficients for the tree 
and penguin operators at the scale mf,. We worked 
with the renormalization scheme independent Wilson 
coefficients. One of the major uncertainties in our 
calculations is due to the fact that hadronic matrix 
elements of both tree and penguin operators involve 
nonperturbative QCD effects. We worked in the fac¬ 
torization approximation, with Nc being treated as 
an effective parameter to include nonfactorizable con¬ 
tributions. We compared our theoretical results with 
the latest experimental data from PDG to determine 
the range of W as 2.07(2.09) < W < 4.54(4.65) 
for q^lml = 0.3(0.5). It has been pointed out that 
the factorization approach is quite reliable in energetic 
weak decays [47] . 

We explicitly showed that the CP violating asymme¬ 
try is very sensitive to W, the CKM matrix elements 
and the form factors. There is a maximum value for 
the differential CP violating parameter when the low 
invariant mass of the tt+tt” pair is near the vicinity 
of the w resonance, 0.780 — 785 GeV. We determined 
the range of the localized integrated CP asymmetry 
with and without p-uj mixing in the ranges of W, the 
CKM matrix elements, and q^/rn^. For all the models 
investigated, we found that the localized integrated 
CP violating asymmetry with p-uj mixing varies from 
-0.0403 (-0.0290) to -0.0752 (-0.0626) corresponding 
to q^/ml = 0.3(0.5) when 0.750 < m{'K'^TT~)iow < 
0.800 GeV. If one adds the events in the cos 9 < 0 
and COS 0 > 0 experimental regions together when 
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0.750 < m(7r+7r“)/ou) < 0.800GeV, the total localized 
integrated CP asymmetry is —0.0294 ± 0.0285. Our re¬ 
sults, especially the signs, agree with the experimental 
data when 0.750 < m(Tr'^Tr~)iow < 0.800GeV. We 
note that the signs are positive in this region when 
p-ui mixing is not considered. This indicates that p-uj 
mixing is vital for A^p to be negative in this region. 
It was shown that the maximum localized integrated 
asymmetry in the range 0.750 < < 0.800 GeV can 

reach -0.0752. We also found that p-uj mixing can make 
the localized integrated CP asymmetries move toward 
the negative direction, and therefore contributes to the 
sign change around the p°(770) peak of m{-K^'K~)iow 
This behavior is independent of the form factor models. 
In our calculations, the difference between the localized 
integrated CP asymmetries in the regions 0.470 < ^/s < 
0.770 GeV and 0.770 < Vs < 0.920 GeV varies from 
0.0076 to 0.0387. Our results by including the p-uj 
mixing mechanism in these two regions and the whole 
region around the p°(770) peak are consistent with 
the experimental results for some models of the form 
factors. 

At this stage, we cannot explain the LHCb exper¬ 
imental data in the regions of positive and negative 
COS0 individually. This is because three-body decays 
of heavy mesons are more complicated than two- 
body decays as they receive more contributions from 
different mechanisms, for example, nonresonants [48], 
the interference between intermediate resonances and 
hnal-state KK i—> tttt rescattering. We will investigate 
the angle distribution of Aqp when considering both 
the p-uj mixing mechanism and the interference between 
different spin intermediate resonances [5]. We will also 
apply more accurate data in the future to further 
decrease the uncertainties in the calculations. With 
parameters with smaller uncertainties, we expect to be 
able to obtain the effects of p-uj mixing more precisely. 
This is important to interpret the angle distribution 
and the sign change of the CP asymmetry around the 
p°(770) peak of m^(7r+7r“)jotu more accurately. 
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